Enhancement of the scissors mode of an expanding Bose-Einstein condensate 
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We study the time-evolution of the scissors mode of a Bose-Einstein condensate during the ballistic 
expansion after release from the magnetic trap. We show that despite the nontrivial character of 
the superfluid expansion, the sinusoidal behavior of the scissor oscillations is recovered after an 
asymptotic expansion, with an enhancement of the final amplitude. We investigate this phenomenon 
with a condensate held in an elongated magnetostatic potential, whose particular shape allows for 
the excitation of the scissors mode. 
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The close connection between Bose-Einstein conden- 
sation in dilute atomic gases and superfluidity has been 
shown by recent research. The most striking signatures 
of superfluidity observed are quantized vortices re- 
duction of dissipative phenomena ^j, scissors modes 
H U , and irrotational flow ||, 0] . 

Due to the particular requirements for the trapping po- 
tential, the observation of the scissors mode has been so 
far reported only for disk shaped condensates [0, where 
the superfluid behavior is evidenced by the different os- 
cillation frequencies with respect to the thermal gas j|. 
Very recently, the authors of || [7j have also analyzed the 
implication of irrotationality on the expansion of Bose- 
Einstein condensates (BECs) initially confined in a ro- 
tating trap, showing that the condensate expands in a 
distinctively different way with respect to a non-rotating 
condensate, thus providing a direct manifestation of the 
quenching of the moment of inertia. 

In this paper we show that the peculiar behavior of 
a free-expanding rotating condensate has also important 
consequences on the evolution of a scissors mode initially 
excited in the trap. Although the condensate under- 
goes a nontrivial expansion, we find that asymptotically 
the scissors mode recovers the initial sinusoidal behavior, 
with an enhancement of the final amplitude. To verify 
our predictions we take advantage of the particular shape 
of the magnetic field in a magnetostatic trap, which per- 
mits the excitations of scissors modes in elongated BECs. 
This allows us to investigate the expansion of a rotating 
superfluid in a novel regime, where the implications of 
the irrotationality are the most dramatic. 

The wave function of a condensate at zero temperature 
can be conveniently written in terms of the density p and 
the phase S 



(1) 



In the Thomas- Fermi approximation the Gross-Pitaevskii 
equation for ip can be transformed in two coupled hydro- 
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dynamics equations for the density p and the velocity 
v = V5/m|| 



dp 
di 



V (pv) = 



(2) 



where U is the trapping potential, g = AnT^a/m the 
coupling strength, m the atomic mass and a the inter- 
atomic scattering length. 

Here we are interested in the solution of Eqs. 
which corresponds to a scissors mode, that is a shape- 
preserving oscillation of the condensate Q. To account 
for this we consider a harmonic potential of the form 
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ij=l 



(t)Xj 



(4) 



where W — W T for symmetry reasons, with initial con- 



ditions Wij(t < 0) 



By tilting the trap potential by a small angle around 
the z-axis at t — 0, one can excite a scissors mode in the 
x — y plane, identified by the angle 



%(t) = oi cos(w sc i), 



whose frequency is u s 



(^+^) 1/2 



(5) 



In this case 



Eqs. (]2|)-(|3|) can be solved exactly with a quadratic ansatz 
for the condensate density and phase |)| |Io| 

p(x,f) = ( Po (t) -\Y, ■i;.\, j :j.,-^J , (6) 

I 1 3 \ 

S(x.,t) = muj ho s Q (t) + - ^ XiBij(t)xj . (7) 
V U =1 / 

Therefore, in general the system can be described by 10 
time dependent dimensionless parameters: po, So, Aij 
and Bij (A and B are 3x3 symmetric matrices, reflecting 
the symmetry property of U, with = = B^ due 



2 



to the vanishing of Wis(t)). These parameters obey the 
following set of first order differential equations [fijl 



dr 
dA 

dB_ 
~d7 



rp O dS 

2ATtB + {A, B} 
W - A-B 2 



(8) 
(9) 
(10) 



with boundary conditions determined by the equilibrium 
distribution of the condensate at r = ujhot — 0. Notice 
that these equations depend on the number of atoms N 
and on the scattering length a only through the initial 
value of p (0) = 0.5a 2 ho (15Na/a ho ) 2 / 5 . 

The subsequent angular motion is evaluated in terms 
of the rotation which diagonalizes the quadratic part of 
the density p(x,t), (that is the matrix Aij). The angle 
of rotation is fixed by the relation M 



tan(20 o ) 



2A, 



An - A 



22 



(11) 



The same Eqs. (§|)-(|h]) can be used to study the ex- 
pansion of the condensate after the release from the trap 
(W = 0), with the appropriate initial conditions. 

Experimentally we study the scissors mode of a BEC 
of 87 Rb atoms in an elongated potential created by a 
QUIC magnetostatic trap |ll]]. The potential has a 
cylindrical symmetry around the x-axis, with frequencies 
w x =2ttx16.3 Hz and ui v =ui z =2ttx200 Hz. The BEC is 
produced by standard radio-frequency evaporation, and 
is typically composed by 10 5 atoms. In all the measure- 
ments we report no thermal component was detectable. 

Since the trap is static, the scissors mode cannot be ex- 
cited with the technique demonstrated for time-orbiting- 
potential traps ||. However, we can take advantage of 
the fact that gravity breaks the symmetry of the mag- 
netic potential resulting in anharmonicity even at the 
trap minimum. Actually, we found that the eigenaxes 
of our magnetic trap are rotated in the x — y plane in a 
localized region around the minimum. The rapid motion 
of the condensate through such deformation is therefore 
equivalent to a sudden rotation of the trap, like that con- 
sidered in Eq. (||). In fact, we found that when we excite a 
dipolar oscillation along the weak a;-axis by a sudden dis- 
placement of the trap minimum, a scissors mode appears 
as soon as the BEC travels for the first time through the 
minimum. Since for an elongated BEC the scissors mode 
proceeds essentially at to y 3> lo x , it is possible to study 
several scissors oscillations before the BEC travels again 
through the center of the trap. By choosing an appro- 
priately small dipolar amplitude, we observe a "pure" 
scissors mode with amplitude #oi — 15 mrad, without 
any apparent shape deformation. 

To describe the characterization of such mode, we need 
first to recall the peculiar behavior of a rotatingBEC dur- 
ing the expansion after release from the trap B. In fact, 
due to the irrotational nature of a superfluid (V x v = 0) 
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FIG. 1: Absorption images of a rotating, elongated Rb con- 
densate during the ballistic expansion. In a) the BEC has 
just been released from the trap (t e x P = 2 ms), and in b) it 
has expanded for 23 ms. The other two pictures, not to scale, 
show the corresponding evolution of a scissors mode. The ini- 
tial oscillation about the horizontal a>axis (c) is transformed 
in an oscillation about the vertical j/-axis for long expansion 
times (d). 



the moment of inertia Q of a condensate is quenched with 



respect to the rigid body value <d r i g |12j 



(x 2 + y 2 ) 



2\2 W «S' 



(12) 



Therefore, due to conservation of energy and angular mo- 
mentum, a rotating condensate cannot reach a symmet- 
ric configuration during the expansion, and undergoes a 
rapid rotation causing the inversion of the "long" and 
"short" axes. In particular, if at release the condensate 
is rotating in the x — y plane, with its long axis forming 
a small angle 9q with the a;— axis, as shown in Fig. [j], ini- 
tially it will start expanding in the short direction. Then, 
when the aspect ratio approaches unity, the condensate 
will continue to expand in the long direction, with its 
long axis rotated by an angle which is close to but not 
exactly tt/2 g). 

We have studied the evolution of the rotation angle 9 of 
the condensate as a function of the expansion time t exp , 
as shown in Fig. |^(a). The measurements are performed 
by releasing the condensate after a fixed evolution time 
t of the scissors oscillations in the trap. For compari- 
son we also plot the angle evolution of the thermal cloud 
just above condensation, which has been put under ro- 
tation with the same technique described above. Note 
that the rotation of the thermal cloud can be experimen- 
tally investigated only until the cloud becomes spherical 
at t exp ~ l/u> x . Initially, the expansion of the BEC and 
the thermal cloud is almost indistinguishable, being char- 
acterized by the angular velocity Qq at release from the 
trap. This is due to the fact that the moment of iner- 
tia of an elongated condensate is close to the rigid body 
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FIG. 2: a) Evolution of the rotation angle for a condensate 
(circles) and a thermal cloud (triangles) released with small 
angular velocity from an elongated trap, b) Evolution of the 
aspect ratio of the rotating condensate. When the aspect 
ratio gets close to unity the angular velocity shows a rapid 
increase. The lines represent the theoretical predictions. 



value (see Eq. (|12|)). The quenching becomes important 
as soon as the aspect ratio approaches unity, where the 
non-classical behavior of the condensate is evident, as 
predicted in Q. Indeed, its rotation undergoes a fast 
acceleration (Fig. ||(b)), followed by a slow evolution to- 
wards an asymptotic angle, close but smaller than ir/2. 
The agreement of experiment and theory is quite good. 
In contrast, the rotation angle of the thermal cloud con- 
tinues to follow the predicted behavior jfj] 



0(t; t exp ) = 6 (t) + ata,n(n Q (t)t exp ) , 



(13) 



as shown in Fig. g(a). This figure shows that the use 
of very elongated condensates allows for the investiga- 
tion of novel regimes where the superfluid expansion is 
dramatically different from the classical one. 

By analyzing the expansion of the condensate for the 
various conditions of initial angle and angular velocity at 
release from the trap, it is possible to reconstruct the evo- 
lution of a scissors oscillation for the expanding BEC. As 
an example, in Fig. we show the calculated evolution of 
the angle during the time of flight, for conditions cor- 
responding to very small angular velocity and maximal 
angular velocity of the trapped oscillations. The behav- 
ior of 6 in all these regimes is qualitatively similar to that 
shown in Fig. ^|. However, in the situations of small ini- 
tial angular velocity the transient rotation is even faster 
than in the other cases, while the asymptotic angle is 
closer to 7r/2. 

By studying the solutions of Eq. (|s|)-(|l0|), we can iden- 
tify three different regimes for the expanding scissors 
mode, depending on the expansion time t exp . i) For 
texp < the condensate expands in a way similar 

to a classical gas, and the sinusoidal behavior of the scis- 
sors mode is preserved, as is possible to see from Eq. (13). 
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FIG. 3: Calculated rotation angle of the condensate during 
the expansion, released at different times during the scissors 
oscillation in the trap, as shown in the inset. Cases A and 
C correspond to small angular velocities, B and D to the 
maximum velocity. The evolution of D is the opposite of B 
(negative angles). 



ii) When t exp ~ 1/uj x the rotation angle is always close 
to ±7r/4, depending on the sign of the angular velocity 
at release. The angle oscillation is therefore close to a 
square wave with the same phase of the angular velocity 
in the trap £lo (t) . iii) For t exp > 1 jio x the angle of the 
condensate approaches 7r/2, and it is more convenient to 
study the oscillations of the angle </> from vertical. We 
find that the original sinusoidal behavior of the oscilla- 
tions is restored, with a time-dependence which can be 
fitted by the form 



(t) + F(t exp )n {t) 



(14) 



= 8 Ql ^Jl+ uj 2 sc F 2 (t exp ) cos(uj sc t + if) , 

where F is a non-trivial function of the expansion time 
and of the trap geometry, and f is a phase shift. This 
behavior is remarkable, since in passing over the "criti- 
cal" time-region at t exp ~ l/u> x , the condensate does not 
lose memory of its initial angular velocity, and it starts 
to behave again as a classical gas. 

We have been able to verify experimentally these ex- 
pectations in our Rb BEC. Since the finite resolution 
of the imaging system prevents us from studying the 
condensate inside the trap, we probe the oscillation of 
the condensate after a minimum expansion t exp = 4 ms, 
as shown in Fig. |4|(a), and we reconstruct the scissors 
mode in the trap. In Fig. ^(b) we compare such a mo- 
tion to the experimental scissors mode after a long ex- 
pansion of 23 ms. The initial scissors mode is clearly 
amplified and phase-shifted after the expansion. For 
the experimental parameters we calculate indeed that 
uj sc F(t exp = 23 ms) ~ 10.7 is much larger than one, 
and therefore the time evolution of the scissors mode 
is substantially a replica of that of the angular veloc- 
ity in the trap. We have observed a similar behavior also 
for collision-induced scissors oscillations of a binary BEC 
system [O. 
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FIG. 4: Evolution of the scissors mode during the ballistic 
expansion. In a) the angle oscillation 8(t) after a short expan- 
sion of 4 ms (circles) is compared to theory (dashed line) to 
extrapolate the oscillation in the trap. In b) we compare the 
latter (continuous line) to the oscillation <j>(t) after a longer 
expansion of 23 ms (circles) . The data show an enhancement 
of the scissors mode, as expected from theory (dashed line). 
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FIG. 5: Calculated amplification factor y/T + ujI c F* of a scis- 
sors mode as a function of U)x/(*>y, for various expansion times: 
texp = 25, 50, 100, 500 ms. Here the the scissor frequency is 
kept fixed to ui ac — 2n x 200.7 Hz and the radial and axial 
frequencies are varied accordingly. 



To study the role played by the trap geometry in 
the amplification mechanism, we have calculated the ra- 
tio between the final amplitude and the trapped one 
as a function of the trap anisotropy u) x /u> v , for vari- 
ous expansion times. In Fig. ^ we show the results, ob- 
tained by keeping fixed the scissor frequency to the value 
2tt x 200.7 Hz, and varying the radial and axial frequen- 
cies accordingly. The behavior shown is independent of 
the initial angle 0Qi, m the small angle regime. This pic- 
ture confirms that in cigar-shaped traps, for typical ex- 



pansion times (t e 



10 -j- 30 ms), the scissors mode 



amplitude is substantially amplified by the expansion. It 
also shows that the amplification factor asymptotically 
tends to a value ~ 2, regardless the trap geometry. It is 
worth noticing that this asymptotic behavior is similar 
to that of a 2D condensate in the limit oj x /u) y — 1, where 
the relation (Q) can be demonstrated analytically |l4) . 
Notice that in contrast Eq. (|l^) would not imply any am- 
plification of the asymptotic oscillations of the thermal 
cloud. 

As a final remark we recall that in general the scissors 
modes are coupled to compressional modes, since a rota- 
tion of the condensate can induce shape deformation if 
the rotation angle is not small enough In our trap, by 
increasing the amplitude of the dipolar oscillation, and 
hence the angular velocity acquired by the condensate, it 
is also possible to access the non-linear regime, where the 
scissors mode couples to quadrupole shape oscillations. 

In conclusion, we have investigated the expansion of 
the scissors mode of BECs. We have shown that the 
asymptotic oscillation is still characterized by a sinu- 
soidal behavior, which depends linearly on the initial am- 
plitude and angular velocity, resulting in an enhancement 
of the final amplitude. The capability of exciting scissors 
oscillations in elongated condensates opens the possibility 
of future studies of rotating superfluids in a novel regime 
of large moments of inertia and small angular velocities. 
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